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FOREWORD 

This  report  presents  an  extension  of  the  ideas  contained  in 
Appendix  III  of  AFFDL-TR-72-83,  Volume  This  effort  extends  the 
formulation  given  in  Appendix  III  tc  compressible  flows  and  accounts  for 
forces  in  non-uniform  flow.  Upwash  and  sidewash  modeling  has  been 
included  for  the  stores  and  ejector  rack  adjacent  to  the  store  for  which 
loading  calculations  are  being  performed.  Comparisons  of  results  with 
experimental  data  are  presented.  This  report  covers  work  performed 
between  February  1979  and  July  1980  by  Calvin  L.  Dyer  and  Charles  B. 
Heath  (AFWAL/FIGC) . 
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SECTION  I 
INTRODUCTION 

This  effort  extends  the  FDL  Subsonic  Store  Separation  Program  by 
changing  the  aerodynamic  models  for  the  separating  and  adjacent  stores  as 
well  as  the  ejector  rack  model.  The  method  presented  here  is  an 
extension  of  the  work  described  in  Appendix  III  of  Reference  1.  The 
advantage  of  this  approach  is  that  it  does  not  employ  the  slender  body 
assumption  of  the  linear  form  of  the  boundary  condition. 
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SECTION  II 

FLOW  TANGENCY  BOUNOARY  CONDITIONS  ON  THE  BODY  SURFACE 

The  coordinate  systems  used  to  calculate  the  source  and  doublet 
strengths  and  associated  velocities,  l.e. ,  velocities  used  In  boundary 
condition  formulation  and  resultant  velocities  once  the  source  and 
doublet  strengths  have  been  determined,  are  shown  In  Figure  1.  The 
transformation  o'  the  non-uniform  flow-field  velocity  vector  from 
cartesian  coordinates  to  cylindrical  coordinates  Is  given  by 

V ( x  ,  r , e )  =  [V(x)cose  +  W{x)sine  +  W(x)sine]  er  +  (1) 

[W(x)cose  -  V(x)sine]  ea  +  U(x)e 

0  X 

where  U(x),  V(x)  and  W(x)  are  the  components  of  the  velocity  vector  In 
the  x,y,  and  z  directions,  respectively. 

If  the  body  surface  Is  specified  by  a  function 

F(x,  r,  e.)  «  0 

then  the  unit  normal  vector  to  the  body  surface  Is  given  by 

>_  _  grad  F  - - - 

n  =  | grad  F| 

Where  grad  F  •  f£  er  +  }  ee  ♦  f|  ex 

The  perturbation  flow  field  Induced  by  the  store  body  can  be  characterized 
by  the  perturbation  velocity  vector  q(x,  r,  e),  i.e,, 

q  (x,r,e)  =  ur(x,r,e)er  +  u e(x,v)ee  +  ux(x,r,e)ex  (4) 

Where  uv,  u^,  and  u.  ar?  the  perturbation  velocities  in  the  x,  r,  and  e 

X  I  o 

directions  (see  Figure  1),  then  the  flow  tangency  boundary  condition  Is 
given  by 

iR,n*o  (S) 


? 


r 
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where  ^  (the  resultant  velocity)  is  the  sum  of  V,  after  the  transformation 
given  by  Equation  1  has  been  applied,  and  q,  given  by  Equation  4. 

Expanding  Equation  5  by  using  Equations  1  and  4  gives 

[  V  ( x )  cos  e  +  W(x)  sin  e  +  u  _  (x.r.s)]  if  + 

r  3r 

[W(x)  cos  e  *  V(x)  sin  e  +  ua  (x,r,e)]V  + 

8  o  Q 

[U (x)  +  u  x ( x  ,  r  ,  e )  ]  |^=0  (6) 

on  F  (  x  ,-r  ,  e  )  -  0  . 


Equation  6  is  the  non-linear  flow  tangency  condition  for  an  arbitrary  body 
immersed  in  a  non-uniform  flow- field.  We  are  going  to  restrict  ourselves 
to  a  closed  body  of  revolution  from  this  point  on  in  this  report.  With 
this  in  mind,  the  local  body  radius  (r)  is  a  function  of  x  alone.  So, 

r  *  R  (  x  ) 

which  can  be  rewritten 

r  -  R  ( x )  =  0 

and  since  F  (x,  r,  e)  =  0  describes  the  body  shape 
r  -  R ( x )  =  F  (x,r,e) 

which  leads  to 

3F  ,  3  F  dR  ( x )  3  F 

3  r  ’  3  x  ~  ‘  d  x  '  ’  3  3  u 

Substituting  into  the  flow  tangency  condition  (Equation  6  yields 

ur  (  x  ,  r ,  e )  =  [  U  ( x )  +  Ux  ( x  ,  r ,  e )  j  fj.(  x. )  _  v  (x)cose-W(x)  s  i  n  e  (8) 

on  the  bcdy  surface. 
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The  boundary  conditions  of  Equation  A  will  be  split  along  the 
cartesian  axes.  This  allows  the  problem  to  be  handled  as  the  super¬ 
position  of  three  separate  problems.  The  three  problems  are: 

(1)  Axlsymmetrlc  flow  past  a  body  of  revolution  with  the  flow  field 
characterized  by  velocity  U  (x).  This  part  of  the  problem  will  be 
modeled  with  three-dimensional  point  sources,  hence  the  subscript  s. 

urs  (  x ,  r )  -  [  l!  (  x )  +  uxs  (x.r)]  (9) 

on  r  =  R  (  x  1  . 


(2)  Lateral  flow  past  the  body  characterized  by  velocity  V(x). 
This  part  of  the  problem  will  be  modeled  by  horizontal  doublets  and  Is 
denoted  by  the  subscript  d,v. 


rd.v  (x'r's!  - 
on  r  = 


xd,v  (x'r’e) 
R  {  x ) . 


dRjxJ 

dx 


-  V  (x)  cos  e 


(3)  Vertical  flow  past  the  body  characterized  by  velocity  W  (x). 
This  part  of  the  problem  will  be  modeled  by  vertical  doublets  and  Is 
denoted  by  the  subscript  D,W. 


rd,w(x,r’e>  = 


'x^w  (x.r.e)  -  W(x)  sin  •  (11) 


jn  r  =  R  (  x  ) 


4 


AFWAL-TR-81 -  3059 


SECTION  III 

TREATMENT  OF  COMPRESSIBILITY  AND  FORMULATION  TO 
DETERMINE  SINGULARITY  STRENGTHS 

At  this  point,  the  present  work  departs  from  that  of  Appendix  III 
of  Reference  1.  The  approach  used  to  solve  for  the  singularity  strength 
is  to  transform  the  problem  to  an  equivalent  incompressible  one,  deter¬ 
mine  the  incompressible  singularity  strengths,  and  then  transform  the 
resultant  flow  field  back  to  the  compressible  case.  In  Reference  2, 
pages  568  through  581,  Karamcheti  discusses  the  various  aspects  of 
determining  point  singularities  for  incompressible  flow  past  a  body  of 
revolution.  The  principal  transformations  used  are 


u'  =p"u,v'-6v  ,  w 1  =  gw 


The  primes  denote  incompressible  space.  The  unprimed  quantities  are  in 
compressible  space  (physical  space).  These  transformations  are  explained 
in  Reference  3,  pages  5  through  7.  By  applying  these  transformations, 
we  can  use  the  various  incompressible  velocity  formulas  given  in  Appendix 
III  of  Reference  1.  In  the  equations  that  follow  in  this  section,  the 
primes  have  been  dropped  for  convencience,  i.e.,  x  in  the  formulas  should 
be  x'  and  r  =  r'  =  (y‘)  +  (z1)  ,  but  since  y'  =  y  and  z'  *  z,  the  prime  on 
r  is  superfluous. 


The  axisymmetric  problem  may  be  modeled  by  a  distribution  of  three- 
dimensional  point  sources  along  the  store  body  axis.  The  incompressible 
velocity  components  produced  by  N  sources,  located  on  the  X  axis  (Figure  1) 
at  x^ ,  1=1,  N,  at  a  field  point  (x,r)  are  given  by: 


(x.r) 


n 


x,)2 


+  r2 


3/2 


x  -  x. 


Xi)2 


+  r 2  3/; 


> 


(14) 
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where  the  Ms^'s  are  the  as  yet  undetermined  source  strengths.  By 
substituting  Equation  14  Into  Equation  9  we  get  N  equations  for  the  N 
source  strengths,  Ms^,  If  we  can  prescribe  uxs  (x,r)  (the  perturbation 
velocity  at  the  N  control  points). 


1 

4  7T 


E 

1=1 


R(x,.)  -  (x.  -  x . ) 


dR(x.) 


dx' 


( 


(xk)  -  X1 )2  +  R2(xk 


m 

) 


3/  : 


=  U(x.) 


dR(xk) 


dx 


(15) 


0  ±  xk  1  ;  k  =  1,2,..  .n-2 


By  using  the  method  described  in  Reference  4,  we  can  define  U(xk,R)  for 
the  aircraft/store  configuration  of  our  choice.  Two  of  the  above 
equations  will  be  deleted  in  order  to  enforce  the  existence  of  a 
stagnation  point  at  the  store  nose,  and  assure  that  conservation  of 
mass  is  satisfied.  By  setting  xk=0,  which  implies  R(xk)-0,  in  Equation  15 
we  force  the  store  nose  to  become  a  stagnation  point. 


u(0,0) 


(16) 


The  conservation  of  mass  is  satisfied 

strengths  equal  to  zero: 
n 


=  0 


by  setting  the  sum  of  the  sources 


(17) 


This  allows  us  to  apply  Equation  15  at  N-2  points  along  the  store  body 
surface.  Experience  has  shown  that  spacing  the  sources  proportioned  to 
the  local  body  radius  and  applying  Equation  15  at  the  midpoint  between 
successive  sources  gives  the  best  results.  If  Equation  15  was  applied  at 
midpoint  between  the  N  sources,  there  would  be  N-l  equations.  Again, 
experience  has  shown  that  if  we  delete  the  fourth  control  point  upstream 
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from  the  point  on  the  body  where  the  maximum  radius  Is  first  reached, 
we  get  a  favorable  solution. 


The  second  or  sldewash  problem  Is  represented  by  a  distribution  of 
three-dimensional  point  doublets  whose  axes  are  aligned  horizontally  to 
oppose  the  sidewash,  i.e.  flow  In  the  y  direction  in  Figure  1.  The 
incompressible  velocity  components  produced  by  N  doublets  located  on  the 
store  bodj  axis  (x  axis  Figure  1)  are  given  by: 


Substituting  Equations  18  into  the  boundary  conditions  given  by 
Elation  10  gives: 


n 


°ixk_<£s;  k  =  1  , 2  ,  .  .  .  n 
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When  xk  denotes  the  control  point  location  and  R(xk)  Is  the  body 
radius  at  xfc.  The  doublets'  locations  and  control  points  are  chosen  to 
be  at  the  midpoint  between  each  source.  This  will  provide  N-l  locations, 
so  the  Nth  doublet  ar.d  control  point  were  arbitrarily  chosen  to  coincide 
with  the  last  source  location.  Reference  4  Is  used  to  determine  V(xk). 

The  third  problem,  upwash.  Is  solved  In  the  same  manner  as  the 
si  dewash  problem.  N  vertical  doublets  are  placed  on  the  store  body 
longitudinal  axis.  The  incompressible  velocity  components  are  given  by: 


u_  ( x  ,  r  ,  e  )  = 
rd ,  w 


sin  e 


n  pd  w 

ip  _ qi  ,w _  _ 

h  [tx  -  «,)*  * 


3  s i n  9  Ip 

4  7T  L~t 


u  .  r 
d  ^ ,  w 

(x  -  Xi)2  +  r2 


u.,  (x.r.e)  = 


u  „  (  x  ,  r  ,  e  ) 

xd  ,w 


cos  e  r 

1  “l  fix  -  x,)!  ♦  r2J  il~‘ 


3  s  i  e  Y"* 

4  n 


"d1.»r,x  '  x<’ 


1-1  [lx  -  Xi)2  +  r^J 


Substituting  Equations  20  into  the  boundary  condition  given  by 
Equation  11  gives: 


n 


0  1  x|(  i  as  ;  k  *  1,2,  ...ii 

The  nomenclature  for  the  vertical  doublets  and  control  point  locations 
is  the  same  as  for  the  horizontal  doublets.  The  method  of  Reference  4  is 
used  to  calculate  W(xk). 
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SECTION  IV 

SELECTION  OF  PRESSURE-VELOCITY  RELATIONSHIP 

Once  the  singularity  strengths  in  Equations  14,  18,  and  20  have, 
been  determined,  they  can  be  used  to  calculate  the  velocity  at  any  point 
on  the  body  surface.  The  surface  velocity  can  be  related  to  the  surface 
pressure  through  various  relationships.  In  Reference  1,  the  incompressible 
Bernoulli  equation  was  used.  This  restr  cted  the  range  of  Mach  numbers 
for  which  their  results  would  be  valid,  but  it  has  the  nice  feature  that 
the  analytical  expressions  for  the  circumferential  pressure  distribution 
could  be  integrated  in  closed  form.  In  Reference  5  and  other  places,  it 
has  been  found  that  the  isentropic  compressible  Bernoulli  equation  relates 
theoretical  prediction  to  experimental  results  better  than  the  theo» etical ly 
consistent  linear  relationship.  In  order  to  extend  the  present  method 
to  compressible  speeds  and  take  advantage  of  this  improved  correlation 
(Reference  5)  the  isentropic  compressible  Bernoulli  equation  was  selected 
from  Reference  6  as  the  pressure-velocity  relationship  for  this  application. 


n 


o  r  v_i  ? 

^  [1  +  -V-  M„  (  1 

yM  V. 

Cl' 


(22) 


The  disadvantage  to  this  approach  is  that  the  circumferential  pres¬ 
sure  distribution  is  not  readily  amenable  to  integration  in  a  closed  form. 
However,  the  pressures  can  be  numerically  integrated  both  circumferentially 
and  longitudinally. 
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SECTION  V 

PROCEDURE  FOR  CALCULATING  BODY  FORCES  AND  MOMENTS 


The  forces  and  moments  on  the  separating  store  are  calculated  by 
integrating  the  pressure  circumferentially  and  longitudinally  on  the 
surface  and  adding  a  term  to  account  for  buoyancy  when  the  store  is  in 
non-uni ^rm  flow.  Aft  of  the  point  of  flow  separation  a  viscous  cross- 
flow  drag  coefficient  is  used  to  calculate  the  forces. 


The  equations  used  to  determine  the  forces  and  moments  acting  on 
the  store  are: 


Cj ,  = 

% 


J*  sep  , 

A  * 


Cp(x,e)  (-sine)  r(x)  de  dX 


o 

2tt 


S6P  f  Mx.e)  (-cose)  r(x)  de  dX 


f  / 

o  0 

X_„  2tt 

r  i 


(xmom  ”  Cn  (-sine)  r(x)  de  dX 


o. 

2  TT 


r  / 


(xmom  "  X)  cp(x»®)  (-cose)  r(x)  do  dX 


>  (23) 


Figure  1  shows  the  sign  convention  for  theta.  The  buoyancy  contribution 
to  the  forces  and  moments  Is  given  by: 
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N(Br)  sR*R 


x)r2 


d(j^) 

dx_ 


From  Reference  1 ,  pages  27  and  28. 


(24  Contd) 


Potential  flow  is  assumed  along  the  body  up  to  a  specified  point. 
After  this  point,  the  forcer,  and  moments  are  calculated  using  crossflow 
drag  theory.  The  equations  for  the  crossflow  drag  forces  and  moments  are: 


(CF) 
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(CF) 


'(CF) 


2Cr 


(CF)  S 


R 

2Cr 


R 

2Cr 


SR  *R 


2Cr 


Vr 
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fvLxh2.  (»1x1)2  w(x)  d x( 
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5ep 
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(  r  dxs 

1  oo  oo  oo 

-*)  x  dxs 


>  (25) 


oo  oo 


,sep 


(X 


mom 


-x)  r  CM  ♦  ^ 


'u  '  u 

OO  QO 


sep 


from  Reference  1,  page  29.  Guidance  is  given  on  page  29  of  Reference  1 
for  determining  the  point  at  which  separation  occurs  for  bodies  in  uniform 
flow.  No  generalized  formula  has  yet  been  devised  to  determine  the  flow 
separation  point  on  bodies  in  non-uniform  flow  fields. 
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SECTION  VI 

PROBLEMS  ENCOUNTERED  AND  RESULTS  OBTAINED  BY  INTRODUCING 
THE  SOURCE-DOUBLET  METHOD  OF  CALCULATING 
BODY  FORCES  AND  MOMENTS 

Method  verification  proceeded  from  the  relatively  simple  case  of 
uni  form  flow  to  the  complex  flow-field  of  a  realis. .c  aircraft/stores 
configuration.  Comparisons  were  made  with  data  from  References  7,  8,  and 
9,  to  demonstrate  the  validity  and  improvement  of  the  basic  approach. 
Flow-field  data  as  well  as  force  and  moment  data  were  used  to  assess  the 
method.  Free-stream  predictions  were  compared  with  data  at  angles  of 
attack  of  2°,  4°,  6°,  and  8°  at  Mach  numbers  of  .25  and  JO.  The 
realistic  aircraft/stores  configuration  was  built  up  from  a  simple  wing- 
body  configuration  by  the  addition  of  various  components,  i,e.  pylon, 
rack,  stores,  and  engine  inlets.  Non-uniform  flow-field  predictions 
were  compared  with  data  at  0°  and  6°  angle  of  attack  at  Mach  number  of  .70. 

Several  check  cases  were  run  for  the  ogive  cylinder  store  (Figure  2a) 
in  a  uniform  flow.  In  implementing  this  procedure,  the  first  problem 
that  occurred  was  negative  pressures.  By  casting  the  boundary  conditions 
for  the  sources  in  terms  of  the  body  slope  (see  Equation  9)  Instead  of 
the  body  radius,  the  zero  streamline  is  not  necessarily  on  the  body 
surface.  The  use  of  body  slope  boundary  conditions  was  found  to  give 
the  "best",  i.e.  smoothest  body  representation,  see  References  1  and  3. 

No  matter  what  is  used  as  the  boundary  conditions,  the  mathematical 
models  (sources  and  sinks)  are  going  to  produce  some  waviness  in  the 
zero  streamline.  The  problem  is  illustrated  in  Figure  2b.  To  resolve 
this  problem  at  each  place  along  the  body  wh'.ie  ,  ;  irface  pressure  was 

required,  the  location  of  the  zero  streamline  was  calculated.  The 
pressure  calculations  were  then  performed  on  the  zero  streamline. 

Once  the  method  was  giving  numerical  results,  several  other  questions 
arose.  First,  how  many  singularities  are  required  to  give  a  meaningful 
answer?  Of  course,  there  is  a  trade-off  between  the  number  of  singular¬ 
ities  and  computational  time.  To  answer  this  question,  runs  were  made 
varying  the  number  of  singularities.  But  before  this  question  could  be 
resolved,  a  problem  with  the  loading  predictions  on  the  aft  end  of  the 
store  was  encountered  as  illustrated  in  Figure  3.  Figure  3  shows  the 
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jp 

distributed  NORMAL  force  -- along  the  store.  Distributed  loads  obtained 
by  integrating  experlmcn tally  measured  pressures  are  indicated  by  the 
points  with  circles  around  them.  The  dashed  line  shows  theoretical 
prediction  based  on  slender  body  theo.y  as  given  by  Equation  46, 

Reference  1.  The  solid  lines  represent  the  theoretical  predictions 
based  on  the  Source-Doublet  method  presented  in  this  report.  Below  the 
horizontal  axis  is  a  silhouette  of  the  store  shape  for  which  the 
calculations  were  performed.  At  the  top  of  ti  ?  figure  are  the  Mach  number 
and  angle  of  attack  for  which  the  calculations  were  presented.  The 
second  line  of  the  heading  gives  the  number  of  circumferential  segments 
the  store  was  broken  into  for  pressure  calculations  and  integration  in 
conjunction  with  the  Source-Doublet  predictions.  Also  on  the  second 
line  is  the  number  of  singularities  used  to  model  the  store,  i.e.  the 
number  of  sources  and  sinks  or  doublets.  The  axial  division  of  the 
body  was  proportional  to  the  body  radius.  This  was  found  to  enhance  the 
fidelity  of  the  volume  modeling  (source-sink)  of  the  body  in  Reference  3. 
The  control  points  for  both  the  source-sinks  and  the  doublets  were 
located  midway  between  the  sources  on  the  body  surface.  The  doublets 
were  located  on  the  body  axis  at  the  same  longitudinal  location  as  the 
control  points.  The  third  line  is  the  number  of  equal  length  body 
segments  the  store  was  divided  into  for  the  slender  body  theory 
calculations.  Just  below  the  header  information  is  the  legend  for  the 
various  symbols  on  the  plot.  To  the  right  of  the  legend  is  the  integrated 
loads  and  measured  loads  from  a  strain-gauge  balance  inside  the  model. 

The  reference  area  is  the  store  maximum  cross-sectional  area.  The 
reference  length  is  the  store  maximum  diameter.  The  moment  reference 
center  was  chosen  to  be  at  the  store  midpoint.  The  pressure  data  was 
obtained  independently  of  the  strain-gauge  balance  data,  so  the  question 
of  how  pressure  lines  bridging  a  balance  affect  the  load  is  not  a  problem 
in  this  data  set.  The  columns  are  normcl  force  coefficient  (Cn)  and 
pitching  moment  coefficient  (r  Note  that  the  last  pressure  data  point, 
i.e.  the  point  at  x  =  .531  ft,  is  on  extrapolated  value.  This  was  done 
so  Simpson's  Rule  could  be  used  to  integrate  the  loadings. 
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Two  problems  are  evident  in  the  aft  end  loading.  The  jump  in 
loading  at  x  =  .43  feet  is  caused  by  shifting  to  the  viscous  cross 
flow  method  of  predicting  forces  as  prescribed  by  Equation  50,  Reference  1. 
Based  on  the  pressure  data,  and  considering  this  same  behavior  was 
observed  at  lower  Mach  numbers  and  angles  of  attack,  it  was  decided  to 
assume  the  flow  was  attached  over  the  whole  length  of  the  body.  This 
assumption  is  used  for  all  subsequent  calculations  in  this  report. 

Although  the  second  problem,  a  large  negative  load,  in  predicting  the 
loadiny  on  the  aft  end  of  the  store  is  evident  in  Figure  3,  it  is  shown 
much  more  graphically  in  Figure  4.  This  problem  results  from  the  way 
the  body  wake  is  being  modeled.  To  satisfy  the  law  of  conservation  of 
mass,  the  body  must  be  closed.  In  References  1,  3,  and  4,  it  had  been 
common  practice  t^  close  the  bodies  by  modeling  the  wake  with  the  same 
shape  as  the  nose.  This  was  our  initial  approach,  but  as  can  be  seen  in 
Figure  4,  the  resultant  velocities  near  the  aft  end  of  the  store  were 
much  higher  than  experimental  results  indicate.  This  problem  was 
corrected  by  extending  the  body  one  diameter  downstream  and  then  closing 
it  with  a  wake  of  the  same  shape  as  the  nose.  Figure  5a  is  typical  of 
the  loadings  produced  by  extending  body.  Based  on  the  results  illustrated 
in  Figure  5a  and  the  other  comparisons  in  the  free-stream  at  other  Mach 
numbers  and  angles  of  attack.  In  all  subsequent  calculations,  the  wake 
is  modeled  by  extending  the  body  one  diameter  downstream  and  closing  the 
body  with  the  same  shape  as  the  nose. 

At  this  point,  we  are  prepared  to  tackle  the  question  of  how  many 
sources  are  required  to  adequately  model  the  loading.  As  an  aside  to 
this  question  an  investigation  into  type  of  integration  and  the  number 
of  circumferential  segments  required  was  also  undertaken.  Simpson's  Rule 
and  Trapezoidal  Rule  integration  were  both  tried.  It  was  found  that 
twenty  circumferential  segments  yielded  converged  solutions  and  that  for 
twenty  circumferential  segments  either  method  of  integration  produced 
essentially  identical  results,  see  Table  1.  The  Trapezoidal  Rule  was 
employed,  since  it  was  easier  to  implement,  faster  and  less  restrictive 
in  point  spacing.  Figure  6  shows  a  plot  of  Cn  versus  the  number  of 
singularities.  It  can  be  seen  that  in  the  range  of  60  to  70  singularities. 
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TABLE  1 

Cn,  Cm,  VERSUS  NUMBER  OF  SINGULARITIES 
Mach  Number  ■  .25  Angle  of  Attack  ■  8.0 
20  Circumferential  Segments 
T  -  Trapezoidal  Rule  Integration 
S  *  Simpson’s  Rule  Integration 


_ - -  — *  -I 

Number  of 

C 

C 

C 

Singularities 

nT 

"t 

ns 

50 

.26133 

.84407 

.26136 

58 

.26563 

.85719 

.26566 

69 

.26721 

.86170 

.26724 

86 

.26640 

.86569 

.26643 

98 

.26659 

.86616 

.26662 

C 

mS 

.84406 

.85718 

.86169 

.86580 

.86615 
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the  solution  becomes  convergent  for  uniform  flow.  All  subsequent  cal¬ 
culations  assume  that  the  60  to  70  singularities  will  yield  converged 
solutions  for  non-uniform  flow-fields  as  well. 

Figure  5  presents  comparisons  for  Mach  numbers  of  .25  and  .70  for 
angle  of  attack  of  2°,  4°,  6°,  and  8°.  These  comparisons  are  for  the 
ogive-cyl inder  store,  Figure  2,  in  a  uniform  flow-field.  The  pressure 
loadings  were  obtained  with  a  single  longitudinal  row  of  pressure  taps 
in  a  meridional  plane.  The  model  was  rotated  through  360  degrees  in 
10  degree  increments.  Then  the  pressures  were  integrated  circumferentially 
to  obtain  the  running  loads  that  are  plotted  in  Figure  5.  The  data  point 
at  x  =  .1386  feet  is  somewhat  questionable.  It  is  clear  from  the  response 
to  angle  of  attack  and  Mach  number  that  the  tap  is  not  plugged.  This 
problem  could  be  caused  by  a  local  imperfection  in  the  body  surface  or 
maybe  a  small  leak  in  the  pressure  line  which  would  be  proportional  to 
the  pressure  at  the  tap.  It  appears  that  the  value  of  DCn/DX  for  this 
tap  should  be  reduced  in  magnitude  to  conform  to  the  trend  of  the  other 
data  points.  It  is  interesting  to  observe  that  reducing  the  value  would 
degrade  the  normal  force  correlation  between  the  pressure  and  strain-gauge 
data  in  every  case  presented  in  Figure  5.  The  reduction  would  improve 
pitching  moment  correlation  between  the  two  types  of  data  for  all  angles 
of  attack  at  Mach  number  of  .25  but  only  at  8  degrees  angle  of  attack  at 
Mach  number  of  .70. 

The  agreement  between  the  two  types  of  experimental  data  and  theory 
is  not  as  good  as  hoped  for.  The  location  of  the  center  of  pressure  was 
calculated  for  both  types  of  data  trying  to  get  some  insight  into  the 
lack  of  data  to  data  correlation.  In  some  instances,  the  variation  in 
center  of  pressure,  based  on  the  data,  was  greater  than  one  body  diameter. 
The  disagreement  between  the  two  types  of  data  is  not  understood. 

Correlation  between  either  theoretical  method  and  either  type  of 
data  was  less  than  excellent.  Both  theoretical  methods  correlate  better 
with  the  loads  obtained  by  integrating  the  pressure  data.  Usually 
within  10  percent  on  normal  force  and  20  percent  on  pitching  moment. 

The  Source- Doublet  method  did  predict  some  lift  carry-over  on  to  the  body 
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aft  of  the  nose  due  to  Its  three-dimensional  nature.  It  appears  that  on 
the  upper  surface  of  the  body  in  the  x  =  .17  to  .20  foot  range  that  there 
may  be  some  recompressions  on  the  upper  surface  which  is  not  predicted  by 
either  theory.  Of  course.  Slender  Body  Theory  could  not  be  expected  to 
predict  this,  since  for  uniform  flow  it  degenerates  into  a  term  which  Is 
only  a  function  of  the  rate  of  change  of  cross-sectional  area.  The 
Source- Doublet  method's  pitching  moment  prediction  correlated  better  with 
the  strain-gauge  data  than  do  the  Slender  Body  Theory  predictions.  Based 
on  this  and  the  fact  that  the  normal  force  prediction  of  the  Source- 
Doublet  method  were  nearly  as  close  to  experimental  values  as  the  Slender 
Body  Theory  it  was  decided  to  continue  on  by  using  the  Source-Doublet 
method  in  conjunction  with  Reference  4  to  predict  loads  in  the  flow- field 
of  a  realistic  aircraft/stores  configurations. 

The  end  product  of  this  investigation  is  to  try  to  improve  store 
trajectory  prediction.  Store  trajectory  prediction  is  a  three-faceted 
problem:  (1)  flow-field  prediction,  (2)  force  and  moment  prediction,  and 
(3)  integration  of  the  equations  of  motion.  The  first  two  points  are 
addressed  in  this  report. 

By  using  the  method  detailed  in  Reference  4  to  prescribe  U(x),  V(x), 
and  W(x)  in  Equations  15,  19,  and  21,  respectively,  the  procedures  out¬ 
lined  in  Sections  III  and  V  can  be  applied  to  predict  loads  in  the  flow- 
fields  of  realistic  aircraft/stores.  The  method  outlined  up  to  this 
point  has  assumed  that  the  flow-field  is  precisely  known,  and  that  the 
only  problem  is  developing  a  method  to  predict  the  forces  acting  on  the 
store.  The  flow-fields  of  real  aircraft/stores  configurations  are  not 
precisely  known  or  predicted.  There  are  many  problems  with  flow-field 
measurement  and  theoretical  predictions.  With  this  in  mind,  both  flow- 
field  and  force  predictions  need  to  be  scrutinized.  Wind-tunnel 
measured  flow-fields  offer  the  principal  yardstick  to  assess  theoretical 
prediction.  Comparison  with  both  flow-field  and  loading  data  will  be 
presented  to  evaluate  the  method  described  in  this  report. 

Figure  7  shows  the  various  model  components  for  which  comparisons 
are  presented.  N1-B2  denotes  the  wing-body  configuration  shown  in 
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Figure  7a.  P(l/3)  stands  for  the  pylon  shown  in  Figure  7b  attached  to 

N1-B2  at  the  1/3  semispan  position.  TER  is  the  triple  ejector  rack  model 
shown  in  Figure  7c.  S2-S3  represent  dummy  store  on  the  shoulder  stations 

of  the  TER.  S2  and  S3  both  have  the  same  shape  as  the  ogive-cylinder 
store  shown  in  Figure  2a.  (All  stores  are  attached  to  the  TER  so  the 
store  midpoint  is  aligned  axially  with  the  rack  centerline.)  The  store 
loading  plots  in  Figures  8,  9,  10,  and  11  have  essentially  the  same 
information  on  them  as  presented  earlier  for  the  uniform  flow-field 
comparisons.  Some  additional  information  is  presented  in  the  form  of 
ticks  on  the  horizontal  axis.  From  left  to  right  they  are:  local  wing 
leading  edge,  rack  nose,  pylon  leading  edge,  pylon  trailing  edge,  aft  end 
of  rack,  and  wing  trailing  edge  These  are  vertical  projections  made 
perpendicular  to  the  fuselage  reference  axis  down  to  the  store  centerline 
in  the  carriage  position.  The  data  presented  in  Figures  8,  9,  10,  and  11 
are  for  a  longitudinal  traverse  along  the  store  centerline  in  the  carriage 
position.  Figures  8  through  11  represent  a  configuration  buildup,  i.e. 
from  a  simplistic  configuration  to  a  realistic  one. 

The  correlation  between  the  integrated  pressure  data  and  the  strain- 
gauge  data  for  the  non-uniform  flow-field  case  presented  are  not  good. 

This  is  consistent  witn  what  was  observed  for  the  uniform  flow-field  cases. 
If  the  pressure  data  can  be  used  as  an  indication  of  what  the  load 
distribution  looks  like,  then  as  would  be  expected,  the  prediction 
deteriorates  as  the  configuration  becomes  more  complex.  It  should  be 
pointed  out  again  that  these  predictions  are  in  the  most  complex  part  of 
the  flow- field,  i.e.  carriage  position.  Dix  in  Reference  10  has  found 
out  that  large  gradients  in  loads  exist  within  a  few  tenths  of  an  inch 
of  carriage  position.  Dix  also  found  that  extreme  fidelity  in  modeling 
the  ejector  rack  hardware  influenced  captive  loads.  The  loading  presented 
in  Figures  8,  9,  10,  and  11  are  dependent  on  the  quality  of  the  flow- field 
predictions  from  Reference  4.  In  Figures  lie,  d,  and  e,  the  U(x)  and  W(x), 
predicted  by  Reference  4  and  non-dimensional i zed  by  the  free-stream 
velocity,  are  compared  to  velocities  measured  in  the  wind  tunnel.  The 
velocities  given  by  the  solid  lines  were  used  to  make  the  predictions 
presented  in  Figures  11a  and  b.  Figures  lie,  d,  and  e  show  that  the 
flow-field  predictions  generally  predict  the  trends  of  the  data  but  the 
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magnitude  is  not  well  predicted.  Flow-field  comparisons  were  made 
corresponding  to  Figures  8,  9,  and  10  also.  They  showed  diminishing 
agreement  between  theory  and  experiment  s  aircraft/stores  configuration 
became  now  complete.  Other  comparisons  simil  to  those  presented  in 
Figures  R,  9,  10,  and  11,  including  flow-fielua,  were  made  and  the  same 
general  conclusion  that  configuration  complexity  causes  agreement  between 
theory  and  experiment  to  diminish.  These  additional  comparisons  will  be 
published  as  a  separate  document.  One  thing  that  needs  to  be  done  that 
was  not  done  is  to  make  predictions  similar  to  those  shown  in  Figures  11a 
through  e  as  a  function  of  vertical  displacement.  These  types  of 
comparisons  would  indicate  where  flow-field  prediction  needs  refinement 
versus  refining  the  method  of  predicting  forces.  Based  on  the  lack  of 
agreement  in  flow-field  prediction  indicated  in  Figures  11c  through  e 
and  supported  by  other  )mparisons,  it  was  decided  to  try  to  improve  the 
flow-field  model  in  Reference  4. 
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SECTION  VII 

FLOW-FIELD  MODELING  EXTENSION 

The  method  of  flow-field  prediction  in  Reference  4  accounts  for 
"primary  interference"  which  is  defined  on  page  2  of  Reference  1. 

Primary  interference  which  can  be  termed  gross  interference  effects 
and  can  be  viewed  as  the  first  term  in  an  iterative  procedure.  The 
primary  interference  approach  fits  somewhere  in  the  spectrum  between 
modeling  each  component,  i.e.  fuselage,  wing,  and  pylon,  as  though  it 
was  in  the  free-stream  and  a  full  blown  panel  or  inter-active  approach 
where  every  component  influences  every  other  component.  This  approach 
was  selected  to  try  to  keep  the  computer  resources  (both  time  and  core) 
down  to  an  economically  manageable  size.  As  part  of  the  primary  inter¬ 
ference  approach,  stores  adjacent  to  the  separating  store  (the  store  on 
which  loads  are  to  be  calculated)  and  the  ejector  rack  are  modeled  by 
sources  and  sinks.  This  models  only  the  volume  effect  and  only  for  the 
free-stream  case.  The  source  and  sink  strengths  are  calculated  for  each 
individual  body  in  a  uniform  free-stream.  Then  these  models  are 
positioned  in  close  proximity  to  represent  the  complete  aircraft/store 
configuration.  This  procedure  is  theoretically  correct  since  the  govern¬ 
ing  equation  is  the  Laplace's  equation  which  is  linear,  so  linear  super¬ 
position  of  solutions  is  perfectly  valid.  The  method  presented  in 
Reference  4  does  not  attempt  to  model  the  upwash  and  si  dewash  caused  by 
adjacent  stores  and  the  ejector  rack.  In  an  attempt  to  improve  both  the 
flow- field  and  force  and  moment  predictions,  a  cross-flow  model  for  the 
shoulder  stores  and  ejector  rack  was  introduced  into  the  method  of 
Reference  4, 

Essentially  the  doublet  equations  described  in  Section  III  were  used 
to  generate  a  model  of  the  cross-flow  produced  by  the  bodies  that  are  in 
close  proximity.  Doublets  were  placed  in  the  horizontal  and  vertical 
direction  along  the  axis  of  each  store  and  the  ejector  rack  similar  to  the 
way  they  were  used  for  the  separating  store  in  Section  III.  The  computer 
program  of  Reference  4  was  modified  to  allow  the  whole  configuration 
(except  the  separation  store  whose  cross-flow  is  being  modeled)  to 
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contribute  to  V(x)  and  W(x)  in  Equations  19  and  21,  respectively.  Once 
the  horizontal  and  vertical  doublet  strengths  were  determined,  then  thei 
contribution  given  by  Equations  18  and  20  could  be  added  to  U(x),  V(x), 
and  W(x)  in  Equations  11,  19,  and  21  to  account  for  the  cross-flow 
induced  on  the  separating  store. 
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SECTION  VIII 

COMPARISONS  USING  THE  EXTENDED  FLOW- FIELD  MODEL 

In  November  1979,  the  Air  Force  Flight  Dynamics  Laboratory  and  the 
Naval  Weapons  Center  sponsored  a  test  at  Arnold  Engineering  Development 
Center  to  study  the  mutual  interference  associated  with  the  carriage  of 
three  MK-83  bombs  on  a  TER  under  the  left  inboard  pylon  on  the  F-4C 
aircraft.  This  wind  tunnel  test  Is  documented  in  References  11  and  12. 
Comparisons  between  theoretical  prediction  and  wind  tunnel  data  of  force 
and  moment  coefficients  are  presented  as  a  function  of  vertical  distance 
from  the  carriage  position.  The  aircraft  configuration  consisted  of  the 
F-4C,  see  Figure  12,  and  the  left  inboard  pylon  (denoted  PI  on  the 
comparison  plots),  see  Figure  13,  and  the  TER  shown  in  Figure  14.  Some 
configurations  will  have  dummy  stores  on  the  shoulder  stations.  The 
dummy  stores  shape  is  labeled  "actual  configuration"  in  Figure  15.  The 
"actual  configuration"  had  to  be  modified  in  the  aft  end  to  accept  the 
balance  and  sting.  The  modified  afterbody  is  also  shown  in  Figure  15. 

The  suffix  MU  after  MK-83  on  the  comparisons  indicates  modified  after- 
body-unfinned,  while  MF  denotes  modified  afterbody-finned.  The  method 
described  in  this  report  only  predicts  forces  acting  on  the  body.  The 
method  described  in  Reference  1  is  used  to  predict  the  fin  forces. 
However,  velocities  induced  by  the  doublets  of  adjacent  stores  and  the 
ejector  rack  were  included  in  the  fin  force  calculations. 

After  the  computer  coding  associated  with  the  flow-field  model 
extension  described  in  Section  VII  had  begun  the  Flight  Control  Division 
of  the  Air  Force  Wright  Aeronautical  Laboratories  decided  to  discontinue 
work  in  the  area  of  store  separation  prediction.  The  comparisons 
presented  here  are  not  to  the  detail  needed  for  thorough  method 
evaluation.  The  results  are  only  to  illustrate  that  the  coding  is 
working,  hopefully  correctly,  and  possibly  provide  the  foundation  for 
someone  else  to  continue  this  investigation. 

All  comparisons  presented  in  this  section  are  for  a  Mach  number  of 
0.60.  One  of  the  main  purposes  of  the  effort  described  in  this  report  is 
to  include  compressibility.  For  a  thorough  method  evaluation  comparisons 
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must  be  made  at  higher  Mach  numbers.  Data  was  taken  at  various  Mach 
numbers  during  the  test  described  In  Reference  12.  The  main  variables  in 
the  comparisons  presented  here  are  parent  aircraft  configuration  and  angle 
of  attack.  The  group  number  on  each  figure  can  be  used  In  conjunction 
with  Reference  12  to  fully  define  all  test  parameters. 

Figures  16,  17,  and  18  present  results  for  a  given  configuration  at 
angles  of  attack  of  1,  5,  and  9  degrees  respectively.  The  configuration 
consists  of  the  F-4C  with  a  TER  rack  mounted  on  the  left  inboard  pylon. 

A  modified  afterbody  MK-83  is  moved  vertically  (perpendicular  to  the 
aircraft  waterlines)  through  the  aircraft  flow-field  from  its  carriage 
position,  i.e.  ZP/D,  to  a  vertical  distance  equal  to  6  diameters 
below  the  carriage  position.  The  first  data  point  is  not  at  exactly 
ZP/D=0.0  because  of  mechanical  interferences  associated  with  the  test 
hardware.  The  configuration  used  for  comparisons  is  the  least  complex 
configuration  tested  for  which  the  method  in  Section  VII  would  be 
applicable,  i.e.  if  the  method  does  not  show  some  improved  prediction 
capability  for  this  case,  provided  the  flow-field  predicted  by 
Reference  4  and  Section  VII  is  reasonable  then  the  method  should  be 
scrapped.  The  legend  for  the  various  symbols  on  the  plots  are  given  in 
the  upper  right-hand  corner  of  each  plot.  The  dots  with  circles  around 
them  are  data.  SBT  stands  for  slender  body  theory  which  denotes  the 
application  of  the  method  of  Reference  4  unmodified.  S.D.  stands  for 
Source-Doublet  and  is  used  to  denote  the  application  of  the  method  outlined 
in  Sections  III  and  IV  in  conjunction  with  the  flow-field  model  of 
Reference  4.  The  suffix  CFLOW  indicates  application  of  the  two  previously 
described  methods  for  force  and  moment  prediction  with  the  addition  of  the 
cross-flow  modeling  described  in  Section  VII  to  the  flow-field  modeling 
of  Reference  4. 

In  all  three  applications.  Figures  16,  17,  and  18,  the  Source-Doublet 
method  Improved  the  correlation  between  predicted  and  measured  moment 
coefficients  (both  pitching  and  yawing  moment)  over  slender  body  theory 
as  described  in  Reference  4  with  or  without  the  flow-field  modeling  of 
Section  VIII.  For  yawing  moment,  the  Source-Doublet  method  with  cross- 
flow  modeling  of  Section  VII  gave  the  best  correlation  between  theory 
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and  experiment  of  all  four  of  the  methods.  The  Source-Doublet  method 
without  the  cross-flow  model  gave  the  best  prediction  for  pitching 
moment.  A  substantial  Improvement  in  both  pitching  and  yawing  moment 
coefficient  prediction  is  obtained  by  applying  the  methods  of  Sections  III, 
IV,  and  VII  in  conjunction  with  Reference  4  over  using  Reference  4  alone. 
Adding  the  cross-flow  model  does  not  necessarily  improve  the  Source- 
Doublet  predictions  as  indicated  in  Figures  16b,  16d,  end  Figure  17b. 

The  Source-Doublet  method  with  and  without  the  cross-flow  model  slightly 
degrades  the  force  coefficient  predictions  shown  in  all  three  figures 
(16,  17,  and  18).  Including  the  cross-flow  model  has  little  effect  on 
the  slender  body  theory  predictions. 

In  Figure  19,  the  cross-flow  model  contribution  to  the  flow-field 
was  terminated  at  ZP/D=2.0.  From  Figures  16,  17,  and  18,  it  can  be  seen 
that  the  cross-flow  contribution  to  the  prediction  is  very  small  below  two 
store  diameters.  Figure  19  is  for  the  most  complex  TER  loading,  i.e. 
two  shoulder  stores  with  fins  are  present.  The  angle  of  attack  of 
parent  aircraft  is  5.0  degrees.  The  Source-Douolet  method  with  cross-flow 
model  does  an  outstanding  job  predicting  yawing  moment,  see  Figure  19d, 
and  does  a  good  job  predicting  normal  force.  Figure  19a.  Pitching  moment 
and  side  force  predictions  are  both  poorly  predicted. 

Figure  20  presents  results  for  a  finned  modified  after  body  store  and 
one  shoulder  store  on  the  ejector  rack.  The  shoulder  store  has  fins  and  is 
on  the  outboard  station  of  the  rack.  The  fins  on  the  shoulder  store  cannot 
be  modeled  by  the  method  described  in  Reference  4.  This  configuration  is  ?.t 
5  degrees  angle  of  attack  and  is  of  intermediate  complexity,  i.e.  it  is  less 
complex  than  the  one  used  to  calculate  the  results  presented  in  Figure  19, 
but  more  complex  than  the  one  used  for  comparisons  In  Figure  18.  As  was 
mentioned  earlier,  the  methods  presented  in  this  report  cannot  predict  fin 
forces.  Therefore,  the  increment  in  force  due  to  the  fins  predicted  by 
Reference  4  has  been  added  to  the  body  forces  predicted  by  the  methods 
described  in  this  report  to  obtain  the  total  forces  acting  on  the  store. 

For  the  force  prediction  all  four  methods  yielded  essentially  the  same 
prediction.  Normal  force  predictions  are  good,  but  side  force  predictions 
leave  a  little  to  be  desired.  The  Slender  Body  and  Source-Doublet  methods 
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both  give  excellent  agreement  with  experiment  when  the  cross-flow  model  is 
included  in  the  flow-field  prediction  method.  The  Source-Doublet  method 
shows  very  good  agreement  with  the  yawing  moment  data  both  with  and  with¬ 
out  the  cross-flow  model.  Slender  Body  theory  gives  very  poor  prediction 
of  the  yawing  moment. 

Figure  21  is  for  the  most  complex  configuration  for  which  data  in  the 
joint  Air  Force/Navy  test.  The  F-4c  is  at  5  degrees  angle  of  attack  with 
two  finned  stores  on  the  shoulder  station  of  the  ejector  rack.  The  Source- 
Doublet  theory  with  the  cross-flow  model  shows  excellent  agreement  with 
the  experimentally  measured  normal  force.  Agreement  within  the  first  two 
diameters  below  the  carriage  positio'  for  pitching  moment  is  poor  for  all 
four  methods  with  Slender  Body  theory  being  slightly  better  than  the 
Source-Doublet  predictions.  Side  force  prediction  is  poor  for  all  four 
methods.  Again  the  Source-Doublet  method  does  a  very  good  job  of  predict¬ 
ing  the  yawing  moment,  certainly  much  better  than  Slender  Body  theory. 

In  all  comparisons,  which  include  several  not  presented,  the  yawing 
moment  prediction  was  consistently  well  predicted  by  the  Source-Doublet 
method . 
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SECTION  IX 
CONCLUSIONS 

The  Source- Doublet  method  presented  here  has  shown  enough  potential 
to  warrant  further  investigation.  The  flow-field  prediction  of  Reference  A, 
with  the  method  described  in  Section  VII,  needs  to  be  studied.  Flow-field 
data  from  che  joint  AFFDL/NWC  test.  References  11  and  12,  could  be  used. 

The  method  seems  to  consistently  predict  moments  better  than  forces. 
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Coordinate  Reference  Frames  Used  in  Source-Doublet  Store  Model 


(a)  Ogive  cylinc_r  store 


(b)  Zero  streamline  superimposed  on 
actual  body  radius. 


Figure  2.  Ogive  Cylinder  Pressure  Store 
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Figure  4.  Short  Wake  Model  Separation  at  1.0  x/£ 


OCN/OX 


AFWAl-TR-81 -3059 


MACH  NO.  =  .2S  ALPHA  =  2.0 

NO.  CIRCUMFERENTIAL  SEGMENTS  20  NO.  OF  SINGULARITIES  65 
NO.  LONGITUDINAL  SEGMENTS  FOR  S.  0.  THEORY  40 


©  ©  © 


PRESSURE  DATA 
SOURCE-DOUBLET 
S-  B.  THEORY 
STRAIN-GAUGE  DATA 


CN 
.071 
.060 
.072 
.  120 


CM 

.  159 
.210 
.239 
.  143 


Figure  5.  Long  Wake  Model  Separation  at  1.0  x/t 
(Contdj  ,  ,  n 

(a)  Alpha*2  Mach  Number*. 25 


31 


DCN/DX 

.00  0.00  1.00  2.00  3.00  4.00 


AFVIAL-TR-81  -3059 


MACH  NO.  =  .25  ALPHA  =  4.0 

NO.  CIRCUMFERENTIAL  SEGMENTS  20  NO.  OF  SINGULARITIES 


LONGITUDINAL  SEGMENTS  FOR 

S.  B.  THEORY 

40 

CN 

CM 

0  PRESSURE  DATA 

.  145 

.385 

SOURCE-DOUBLET 

.135 

.436 

S.  B.  THEORY 

.  144 

.478 

STRAIN-GAUGE  DATA 

.169 

.379 

i 


-,-j - - 1 - 1 - I  I  * 

*0.00  0.08  0.17  0. 25  0.33  0.42  0-50 

X. FEET 


Figure  5  (Contd).  Long  Wake  Model  Separation  at  1.0  x/s, 
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Figure  5  (Contd).  Long  Wake  Model  Separation  at  1.0  x/t 
(d)  Alpha=8°  Mach  Number=.25 
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MACH  NO.  =  .70  ALPHA  =  2.0 
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Figure  5  (Contd).  Long  Wake  Model  Separation  at  1.0  x/Jt 
(e)  Alpha-20  Mach  Number=.70 
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MACH  NO.  =  .70  ALPHA  =  4.0 
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Figure  5  (Contd).  Long  Wake  Model  Separation  at  1.0  x/d 
(f)  Alpha=4°  Mach  Number=.70 
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Figure  5  (Contd).  Long  Wake  Model  Separation  at  1.0  x/s. 
(g)  Alpha=6°  Mach  Number=.70 
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MACH  NO.  =  .70  ALPHA  =8.0 

NO.  CIRCUMFERENTIAL  SEG  NTS  20  NO.  OF  SINOULARl TIES  65 
NO.  L0NGITU01NAL  SEGMENTS  FOR  S.  8.  THEORY  40 


CN 

CM 

PRESSURE  DRTfl 

.306 

.772 

S0URCE-00U8LET 

.261 

.852 

S.  8.  THEORY 

.287 

.953 

STRAIN-GAUGE  OATA 

.338 

.767 

o 


Figure  5  (Concluded),  long  Wake  Model  Separation  at  1.0  x/t 
(h)  Alpha=1°  Mach  Number=.70 
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Figure  6.  C  as  a  Function  of  Number  of  Singularities 
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Figure  7.  Wind-Tunnel  Models  Used  in  Sample  Calculations 
(a)  Wing-Fuselage  Combination 
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PARENT  CONFIGURATION  N1-B2  Z1=ZER0 
MACH  NO,  =  .70  ALPHA  =  6.0 
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Figure  8.  Wina-Body  Forces  and  Moments  (a)  ^n/dx  vs  x 
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Figure  9.  Wing-Body-Pylon  Forces  and  Moments 
(a)  dCn/dx  vs  x 
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PARENT  CONFIGURATION  N1-B2-P11/3)  Z*i=ZERO 
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Figure  9.  Wing-Body-Pylon  Forces  and  Moments 
(Concluded)  (b)  dCy/flx  x 
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PARENT  CONFIGURATION  N1-02-PI  1/3  )-*TER  Z1=ZER0 
MACH  NO.  =  .70  ALPHA  =  6.0 
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figure  10.  Wlng-Body-Pylon-TER  Forces  and  Moments 

(a)  dCn/dx  vs  x 
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PARENT  CONF  IGURAT  ION  N1-82-PI  1/3  l-TER  Z1--ZER0 
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rigure  10.  Wing-Body-Pylon-TER  Forces  and  Moments 
Concluded) 

(b)  dCy/dx  vs  x 
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PARENT  CONFIGURATION  N1-B2-PI 1/3 J-TER-S2-S3  Z1*ZER0 
MACH  NO.  =  .70  ALPHA  =  6.0 
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Figure  11.  Wing-flody-Pylon-TER-Shoul  der  Stores  Forces  and  Moments  and 
Flow-Field  Data  (a)  c*('n/dx  vs  x 
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PARENT  CONF IOURAT ION  Nl-BZ-Pl 1/3 1-TER-S2-S3  Z1=ZER0 
MACH  NO.  =  .70  ALPHA  =  6.0 
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Fifcontd)  '  w^n9'Boc*/"p/^on'T^R'^^ou^er  Stores  Forces  and  Moments  and 
Flow-Field  Data  (c)  u  Flow-Field  Velocity 
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PARENT  CONFIGURATION  N1-B2-PI 1/3  J-TER-S2-S3  Z1=ZER0 
MACH  NO.  =  .70  ALPHA  =  6.0 
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Figure  11.  Wing-Body-Pylon-TER-Shoulder  Stores  Forces  and  Moments  and 
(Contd) 

Flow-Field  Data  (d)  v  Flow-Field  Velocity 
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Figure  11.  Wing-Body-Pylon-TER-Shoulder  Stores  Forces  and  Moments  and 
(Concluded) 

Flow-Field  Data  (e)  w  Flow-Field  Velocity 
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Figure  12.  F4-C 
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Figure  13.  F4-C  Pylon  Models 
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Figure  19.  F-4C-Inboard  Pylon-TER-2  Finned  MK-83  Shoulder  Stores 

Loads  at  Alpha  .  =4°  with  the  Modified  Unfinned  MK-83 

K  store 

(a)  C  vs  Zn  ,n 


C 


.40  -1.20  -1.C0  -0.00  -0.6C 


AFVIAL-TR-81  -3059 


GROUP  241 
„  .  .OflTfl 


,S .  0 . 

tSBT  CFLOW 
,S .  D .  CFLOW 


3.00 

ZP/O 


5.00  6.00 


F-4C . PI .T. (S2)F, ( S3  )F 


MK-83 . MU 
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Figure  21.  F-4C-Inboard  Pylon-TER-2  Finned  MK-83  Shoulder  Stores 
Loads  at  Phastore=4°  with  Modified,  Finned  MK-B3 
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